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Magnet oelectric MnPSa as a new candidate for ferrotoroidicity 
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We have revisited the magnetic structure of manganese phosphorus trisulfide MnPSs using neutron 
diffraction and polarimetry. MnPSs undergoes a transition towards a colhnear antiferromagnetic 
order at 78 K. The resulting magnetic point group breaks both the time reversal and the space in- 
version thus allowing a linear magnetoelectric coupling. Neutron polarimetry was subsequently used 
to prove that this coupling provides a way to manipulate the antiferromagnetic domains simply by 
cooling the sample under crossed magnetic and electrical fields, in agreement with the non-diagonal 
form of the magnetoelectric tensor. In addition, this tensor has in principle an antisymmetric part 
that results in a toroidic moment and provides with a pure ferrotoroidic compound. 



Recent years have seen a renewed interest for mag- 
netoelectricity and a fast increasing number of stud- 
ies devoted to multiferroics materials [l|. These com- 
pounds are defined as presenting at least two ferroic 
orders amongst the three widely known: ferroelastic- 
ity, ferroelectricity and ferromagnetism. Although their 
antiferro- counterparts seem to be less interesting for pos- 
sible technological applications, they are also often con- 
sidered in the above classification. Recently, a fourth 
ferroic order, ferrotoroidicity, has been proposed as an 
order parameter associated with the breaking of both 
time reversal and space inversion whereas ferroelastic- 
ity is symmetric with respect to these two operations, 
ferromagnetism only breaks the time reversal and ferro- 
electricity the space inversion jlM6|. 

Ferrotoroidicity is intrinsically associated with the lin- 
ear magnetoelectric (ME) coupling from which a mag- 
netization (electric polarization) can be proportionally 
induced by an electric (magnetic) field [7[. The toroidic 
moment can be seen as generated by a vortex of magnetic 
moments and is defined as \^^ri x Mi (with i labelling 
the magnetic atoms in the unit cell). It is proportional 
to the antisymmetric part of the ME tensor a^^, hence 
present only when this tensor has non-diagonal terms ^] . 

Multiferroics usually exhibit complex magnetic orders, 
often driven by magnetic frustration, with e.g. cycloids 
or incommensurate modulated phases. However, for a 
linear ME coupling to be allowed in a centrosymmet- 
ric crystal, a necessary condition is a magnetic struc- 
ture with a zero propagation vector which loses the in- 
version center below the ordering temperature [7[. A 
number of oxides have been identified to fulfill this con- 
dition and to be potential magnetoelectrics [9] . The most 
famous and studied example is Cr203, with a linear diag- 
onal ME tensor associated with a collinear antiferromag- 
netic order without any spontaneous electric polariza- 
tion 10Ml2J. The transition temperature is higher than 
room temperature which is promising for potential ap- 
plications: for instance, the switching of the exchange 



bias has been proven in a Cr203/(CoPt)3 antiferromag- 
netic/ferromagnetic heterostructure. Performing ME an- 
nealing [12], i.e. cooling through the Neel temperature 
in electric {E) and magnetic {H) fields (parallel or anti- 
parallel), the antiferromagnetic domains population was 
unbalanced in the magnetoelectric layer [13]. In Cr203, 
the antiferromagnetic domains are the two so-called 180° 
domains with moments in opposite directions (time re- 
versal conjugates). They have different energies when 
submitted to an electric and magnetic field, the differ- 
ence being proportional to the corresponding elements of 
the magnetoelectric tensor AG = 2a^^ EjHj. The work 



of Borisov et al. on exchange bias pj underlines the im- 
portance of domains which form upon a phase transition 
to an (anti)ferroic phase, and of their manipulation. 

Up to now, two techniques have proven to be able to 
probe such 180° antiferromagnetic domains: second har- 
monic generation optical spectroscopy [145 ^^^ neutron 
polarimetry [15|. This last technique measures the do- 
main population and has been used for instance in Cr203, 
MnGe03 and LiCoP04 to evidence the imbalance of do- 
mains through cooling in E and H 15|, ll6| . 

In the search for new magnetoelectric candidates, we 
have investigated MnPS3 by neutron diffraction and neu- 
tron polarimetry. Its main structural and magnetic prop- 
erties are known since the 1980s |17H19|. but this com- 
pound was never reported to show any magnetoelec- 
tric coupling. MnPS3 (manganese phosphorus trisul- 
fide, sometimes referred to as manganese thiophosphate) 
belongs to a family of compounds mainly studied for 
their intercalation properties [20]. It is a lamellar com- 
pound, crystallizing in the monoclinic space group C2/m. 
MnPS3 is a highly resistive broad band semiconductor 
with a gap close to 3 eV and is optically transparent 
with a green colour [21]. The transition metal ions Mn^+ 
with spin 5/2, responsible for the magnetic properties, 
form a honeycomb lattice in the (a, 6) plane. The weak 
inter layer coupling mediated by the S atoms, was pro- 
posed to be purely of Van der Waals origin. However, the 



analysis of the transition occurring at 78 K towards an 
antiferromagnetic collinear phase, suggests an interplane 
exchange associated with some degree of metal- hgand co- 
valency 0,l23|. 

In this article, we report on the evidence and charac- 
terization by neutron polarimetry of linear ME coupling 
in MnPSs and of its influence on the antiferromagnetic 
domain populations. We discuss these results, in partic- 
ular the consequence of the symmetry of the ME tensor 
with respect to ferrotoroidicity. 

A MnPSs crystal, grown by the method reported in 
reference [22| , has been studied by neutron diffraction at 
the Institut Laue-Langevin. The crystal and magnetic 
structures were first checked on the four-circle CEA-CRG 
diffractometer D15 using a wavelength of 1.174 A, at low 
temperature, in a closed-cycle refrigerator. The crys- 
tal structure was confirmed to be C2/m with unique b 
axis with cell parameters a=6.05(l) A, b=10.52(3) A, 
c=6.80(2) A, a=7=90 ° and /3=107.3(2) ° at 90 K. The 
extinction turned out to be completely negligible. 

Spherical neutron polarimetry experiments were per- 
formed on the hot neutron beam diffractometer D3 with 
a wavelength of 0.84 A, using CRYOPAD. The technique 
consists in setting the polarization of the incident beam 
along different directions, and to analyze the polarization 
vector of the scattered beam. More detailed descriptions 
of the technique can be found in ref. [15|, |l6| . Such mea- 
surements usually yields for each Bragg reflection a 3x3 
matrix P^j, with i, j = X, F, Z, that relates the scattered 
polarization to the incoming one. For each reflection, 
X is parallel to the scattering vector, Z is vertical and 
Y completes the right handed set. The Pij matrix el- 
ements contain all the information on the arrangement 
of the magnetic moments and on the different magnetic 
domains if present. For this experiment the crystal was 
mounted with the a axis vertical (defining Z), to perform 
measurements in the (6*, c*) scattering plane. 

In a first step, the magnetic structure was refined from 
the integrated intensities collected on D15 at 10 K. The 
resulting magnetic structure was found to agree with the 
published one [iSj. Below Tat, the antiferromagnetic 
phase is characterised by a zero propagation vector. The 
four magnetic atoms (0, 0.3327, 0), (0, 0.6673, 0), (0.5, 
0.8327, 0) and (0.5, 0.1673, 0) on the 4g Wyckoff site have 
their magnetic moment coupled following a H \ — se- 
quence, i.e. each Mn^+ is antiferromagnetically coupled 
with its nearest neighbors in the (a, b) plane, and the 
coupling between adjacent planes is ferromagnetic. The 
magnetic moments were reported to lie along c* [18|. 

This magnetic structure description was further sub- 
mitted to the stringent test of CRYOPAD. The polariza- 
tion matrices measured for a set of (0,/c,/) reflections at 
2 K without any ME annealing were refined. For each 
reflection, the 9 Pij and the 9 P-ij terms (obtained by 
flipping the incident polarization) were recorded during 
the experiment, thus eliminating possible systematic er- 
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FIG. 1: (color online) Observed versus calculated polarization 
matrix elements for the published magnetic structure (left) 
where the magnetic moments are along c*, and the revisited 
one (right) where the Mn^^ moments lie at a finite angle 
of ^ —8° from this c* axis. The corresponding magnetic 
arrangements are shown in the bottom panels. 



rors. The resulting 99 observations were fitted using the 
magnetic structure described before. The domain popu- 
lations and the Mn^+ magnetic moment amplitude were 
the only parameters of the fit. The result is shown in the 
left column of Fig. [H In a second step, the angle be- 
tween the moment direction and the c axis in the {c^a'') 
plane was let free to vary. After refinement, was found 
equal to -25.4(2)°: taking into account the P monoclinic 
angle, the moments form an angle of ~ —8° with respect 
to c* instead of being parallel to it. The result is shown 
in the right column of Fig. [H This small change drasti- 
cally improves the fit. The resulting magnetic structure 
is shown in Fig. [21 The magnetic moment amplitudes 
are 4.43(3) /iB, slightly smaller than the expected satu- 
rated value of 5 /i^, but in good agreement with previous 
neutron works p^]. In such a structure, there are only 
two 180° antiferromagnetic domains. Their volume frac- 
tion are vi = 55.3(3)% and V2 = 44.7(3)%, i.e. close 
to an equirepartition. Usually, the domain repartition 
is quantified through the normalized volume difference 



^1+^2 



0.106(4) [15]. The existence of a tilt of the 
magnetic moments with respect to natural symmetric di- 
rections (crystal axes) indicates the presence of compet- 
ing anisotropics, result also suggested by a gap observed 
in the spin waves dispersion curves [22|. It could be due 
to dipolar and single- ion anisotropics, as found in other 
systems with Mn^+ ions [24]. 

The MnPSa magnetic point group is thus 2ym where 
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FIG. 2: (color online) Antiferromagnetic structure on the hon- 
eycomb lattice of Mn^^ . 



the inversion center is combined with time reversal. It is 
of particular interest because it allows a linear ME effect 
J9[. Moreover, the ME tensor is constrained by symmetry 
to be non diagonal with non-equal off-diagonal elements. 
In the (a, 6, c) frame, this tensor is: 
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Given the possibility of a linear ME effect, we tried 
to characterize it, through its influence on the domain 
populations, by neutron spherical polarimetry. Indeed, 
in such a magnetic structure, the nuclear and magnetic 
structure factors are in phase quadrature (real structural 
part, imaginary magnetic part). A spherical polarization 
analysis using CRYOPAD is able to give information on 
the respective proportion of the two antiferromagnetic 
domains. The complete formalism is described in ref. 
[l5|- Calculations show that the matrix elements Pxz^ 
Pzx^ PxY^ and Pyx are highly sensitive to the domain 
population parameter r]. 

In this experiment, the MnPSa crystal has been sub- 
mitted to electric and magnetic fields. The electric field, 
with voltages up to 2 kV, was obtained via two Al elec- 
trodes separated from 7 mm surrounding the sample 
within the cryostat. It could only be applied along c* 
due to the lamellar sample shape. The tail of the cryo- 
stat could be placed within the gap of an electromagnet 
delivering magnetic fields as high as 1 T, parallel or per- 
pendicular to the E field direction i.e. along c* or 6* 
respectively (see Fig. [3|). For these two relative orien- 
tations of E and H., the following combinations of ME 
tensor elements have been probed: 
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TABLE I: Relative domain population 77, after ME annealing 
for different orientations and amplitudes of the electric and 
magnetic fields. The first two lines corresponds to no ME 
annealing and annealing in E only. The orientations of the 
fields are given by the unit vectors c* and h* parallel to the 
c* and 6* axes. 
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with ^ = /3— ^ = 17.3°. Several ME annealings were pro- 
duced by applying different electric and magnetic fields 
on the sample during cooling from the paramagnetic state 
(90 K) to the ordered phase (50 K). The fields were 
then switched off and the cryostat was positioned within 
CRYOPAD for the neutron polarization measurements. 
This procedure was repeated for each (^, H) set. The 
first measurements performed on several {O^k^l) reflec- 
tions without ME annealing allowed to calculate precisely 
the domain population and to fix the different parame- 
ters. For subsequent measurements, only the Pxz matrix 
element of the (0,2,0) reflection was measured as the most 
sensitive term to domain populations (see Fig. [3]). The 
domain populations were determined for increasing val- 
ues of the product EH^ for the parallel or perpendicular 
relative orientations of E and H. 

It turned out that the domain distribution could be 
modified from ~50% of each domain without ME anneal- 
ing to an almost single domain sample under the simul- 
taneous action of E and H (results summarized in table 
[Hand Fig. [3j). In the perpendicular geometry, the sample 
becomes rapidly single domain for fields product equal 
to -0.025 Jm-2ps-^ (typically 0.35 T and 86 kV/m). A 
single domain was also achieved for the parallel geome- 
try but for field products ten times larger (typically 1 T 
and 270 kV/m). According to the symmetry of the ME 
tensor, only a^^ should be non-zero, and therefore act- 
ing on the domain population. The observed imbalance 
of the domain populations in the parallel geometry can 
^^actually be explained by a misalignment of both fields of 
^^ the order of 5°. The apparent dispersion of points for 
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FIG. 3: Upper: Absolute value of the relative domain pop- 
ulations ?7 as a function of the product of the electric and 
magnetic fields for different respective orientations. Lower: 
polarization matrix elements of the (0 2 0) reflection for equi- 
populated antiferromagnetic domains (red dots) and for al- 
most single domain 1 (blue triangles) or single domain 2 
(green squares). 



E J_ H in Fig. [3] is of the same origin. Note that the 
reversal of either E or H favors a selection of the other 
antiferromagnetic domain. Finally, it has been checked 
that no ME annealing was obtained with an electric field 
alone. This experiment thus reveals the particular form 
of the linear ME tensor of MnPSs, in agreement with the 
symmetry of the magnetic space group. 

It is worth noting that this non-diagonal ME ten- 
sor allows in principle ferrotoroidicity, i.e. the align- 
ment of toroidic moments. In the experiment geometry, 
it was impossible to investigate the antisymmetric part 
of the ME tensor proportional to the toroidic moment 



T/e oc {afj — a^f )eijk with Cijk the Levi-Civita sym- 
bol, as it was directly probed in GaFeOa [25]. However, 
a significant volumic toroidic moment could be calcu- 
lated in the monoclinic unit cell: T=(0.0175, 0, -0.0027) 
/i^/A^ [5]. MnPSa thus provides with the first example 
of pure ferrotoroidicity (not associated with ferromag- 
netism and/or ferroelectricity) [26] with intrinsic parity- 
time reversal symmetry breaking. Ferrotoroidic domains 
should be present in MnPSs coinciding with the antifer- 
romagnetic ones and manipulated conjointly [6]. 

To summarize, a non-diagonal linear magnetoelectric 



effect has been established in the quasi 2D antiferromag- 
net manganese phosphorus trisulfide using polarized neu- 
trons with spherical polarization analysis. The associ- 
ated behavior is in agreement with the symmetries of 
the magnetic point group and also compatible with fer- 
rotoroidicity. The ME domains population can be easily 
manipulated under reasonably weak fields through ME 
annealing. MnPSs, with its simple magnetic structure 
and rather high ordering transition temperature, offers 
then a new model system to better understand magneto- 
electricity and ferrotoroidicity. 

We warmly thank A. Cano, E. Katz and P.J. Brown 
for fruitful discussions, the latter also for the help with 
data analysis, and S. Vial for the technical support. 
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